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Abstract. We review an approach for the computation of Feynman integrals 
by use of multiple polylogarithms, with an emphasis on the related criterion of 
linear reducibility of the graph. We show that the set of graphs which satisfies 
the linear reducibility with respect to both Symanzik polynomials is closed 
under taking minors. As a step towards a classification of Feynman integrals, 
we discuss the concept of critical minors and exhibit an example at three loops 
with four on-shell legs. 



1. Introduction 

In recent years we witnessed rapid progress in tlie developement of techniques 
for the computation of higher order corrections in perturbative quantum field the- 
ory. While other talks at this conference cover progress in the computation of entire 
amplitudes, our talk refers to the 'classical' approach of computing the amplitude 
by its Feynman graphs, which is inevitable when meeting the needs of present col- 
lider experiments. In this field of research, it has shown to be fruitful to discuss 
Feynman integrals in their own right, without restrictions to a particular quantum 
field theory. 

Computations of higher order corrections to observables often start from the 
consideration of hundreds or thousands of Feynman integrals with tensor structure, 
and proceed via effective standard procedures to reduce the problem, possibly to a 
relatively small number of scalar integrals. At higher loop-orders, the evaluation of 
the latter remains to be the hard part of the problem. There is no algorithm which 
would succeed in the analytical computation of every Feynman integral. However, 
there is a variety of powerful methods which have been useful for a wide range 
of relevant cases, such as the Mellin-Barnes approach (see [9l I60L 1591 [56]), the 
expansion of hypcrgcomctric functions |48l I39L I40j . differential equation meth- 
ods [43l [53l [33, 49, 50 , difference equations [Ml [57l [58l [13 or position-space 
methods |28| (also see |37| ). In this talk we focus on the approach of iteratively 
integrating out Feynman parameters by use of multiple polylogarithms. 
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In order to choose an appropriate strategy for the computation of a given 
Feynman integral, it would be desirable in general, to know in advance, which are 
the classes of functions and numbers the integral may evaluate to. As a slightly 
more refined question of this type we may ask: Which scalar Feynman integrals can 
be expressed by multiple polylogarithms and multiple zeta values, and for which 
integrals do we need a wider range of functions and numbers? In the past few years, 
questions of this type turned out to define a fruitful common field of research for 
quantum field theorists and algebraic geometers alike. While the physicist's interest 
in these questions is given by the desire to compute specific integrals or to learn 
about the 'number content' of a given quantum field theory, the mathematician 
arrives at the same question from a different direction. In a very general context, 
Feynman integrals can be viewed as period integrals, and the question of evaluating 
to multiple zeta values is related to the question wether an underlying motive is 
mixed Tate over Z (see [121 [2011211 l24l [3ll4] ) . 

A definite classification of Feynman graphs with respect to the above questions 
is missing. However, for vacuum and two-point graphs important progress was made 
by considering the first Symanzik polynomial, given by the Feynman parametric 
representation of a Feynman integral. Even though many vacuum-type Feynman 
integrals evaluate to multiple zeta values |17l I18| . this is not the case in general. A 
first vacuum graph whose period has to belong to a set of numbers beyond multiple 
zeta values was exhibited in a recent article by Brown and Schnetz '24j (also see 
[22j ). When allowing the Feynman integrals to depend on kinematical invariants 
and particle masses, we can ask for graphs where multiple polylogarithms are not 
sufficient to express the result. Here the first cases show up at much lower loop- 
order, such as in the case of massive sunrise graphs and related graphs with a cut 
through three massive edges (see e.g. [6], [7])- 

In this talk we review a criterion on graphs which is related to the above 
questions and show that if a graph satisfies the criterion, its minors do so as well. 
In graph theory such a minor monotony is an important and desireable feature. In 
section 2 we begin with a brief reminder on scalar Feynman integrals, their two 
Symanzik polynomials and the approach of integrating out Feynman parameters 
by use of multiple polylogarithms. In section 3 we briefly review the criterion of 
linear reducibility of a graph, which is used to decide whether a given integral 
can be computed by use of the method. If this is the case, the functions and 
numbers in all intermediate steps and in the result will not exceed combinations 
of multiple polylogarithms and their values at rational points. In this way the 
criterion and the corresponding algorithm are useful tools for adressing the above 
questions. In the case of integrals only involving the first Symanzik polynomial, 
the criterion was extensively studied in [20|, 121] . As the iterated integration over 
Feynman parameters can be expected to be useful in the case of integrals depending 
on kinematical invariants and particle masses as well, we intend to extend the 
discussion to the second Symanzik polynomial. In section 4 we consider linear 
reducibility with respect to both Symanzik polynomials and show that the set of 
linearly reducible graphs is closed under taking minors. This property is useful for 
a classification, as it allows us to characterize families of reducible graphs by a small 
number of graphs not belonging to the family. In a case study we exhibit such a 
'forbidden minor' at the level of massless three-loop graphs with four on-shell legs. 
Section 5 contains our conclusions. 
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2. Multiple Polylogarithms and Feynman Integrals 

In this section we recall some general facts about Feynman integrals, Symanzik 
polynomials and a method to compute period integrals by use of multiple polylog- 
arithms. Let us begin with a generic Feynman graph G with n edges, loop-number 
(i.e. first Betti number) L>1 and with r external half-edges (or 'legs'). We label 
each edge by an integration variable ai (Feynman parameter), an integer Ui (ex- 
ponent of the Feynman propagator), a real or complex variable (particle mass). 
Each leg is labelled by a vector pj (external momentum) . 

To this labelled graph G we associate the scalar Feynman integral in Dimen- 
sional Regularization: 

where u = "Y^^^iVi- (We omit to write a trivial prefactor by which the integral 
becomes independent of the physical mass-scale.) The Feynman integral Iq and the 
function J-q depend on the particle masses and on certain kinematical invariants, 
which are quadratic functions of the external momenta. The functions Ug and J-q 
are the first and second Symanzik polynomial of the graph. A definition is given 
below. Usually a Feynman integral is associated to a Feynman graph by Feynman 
rules in momentum or position space, and we refer to the literature |41l 151] for the 
standard computation leading from there to the Feynman parametric representation 
given in eq. 12.11 

Eq. [211] defines a very general class of integrals which deserves our attention for 
several reasons. Firstly, the class contains the Feynman integrals of scalar quantum 
field theory such as (f)^— or 0^— theory. Secondly, any Feynman integral with a 
tensor-structure, arising from a physical quantum field theory, can in principle be 
expressed in terms of scalar integrals of the above class |57l 158] . Thirdly, as we 
allow the Vj to take arbitrary integer values, there are well-known identities between 
these scalar integrals which can be used for efficient reduction procedures |29| . As a 
consequence, integrals of the above class appear in a wide range of physical set-ups 
and their evaluation is the bottleneck of many computational problems in particle 
physics. 

The parameter D can either be fixed to the integer space-time dimension or, as 
the integral is very often ill-defined in the desired dimension, one may consider Ig in 
Dimensional Regularisation where _D is a complex variable. Then, in order to sep- 
arate the pole-terms and obtain finite contributions in four-dimenional Minkowski 
space, one usually attempts to compute the coefficients of a Laurent-expansion 

oo 
j=jo 

with Z) = 4 — 2e, to a desired order. Even though the computation of the functions 
Cj can be very difficult, we can make a general statement about them. It is shown in 
|15j that if for an arbitrary Feynman graph we evaluate any function Cj at algebraic 
values of the squared particle masses mf and kinematical invariants Si, where all 
mf > and all Si < 0, we obtain a period according to the definition of Kontsevich 
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and Zagier |42j . For the special case where the Feynman integral takes the form 



this statement was already proven in [8 . It seems that Feynman integrals in fact 
evaluate to a restricted subset of periods and it is an important challenge to un- 
derstand which one this is. 

Let us now recall the definition of the Symanzik polynomials of a Feynman 
graph G. The first Symanzik polynomial is defined as 

= X! n 

where the sum is over all spanning trees of the graph G. The second Symanzik 
polynomial is defined as 



i=l 



with 




Here the sum runs through all spanning two- forests (Ti, T2) of G, where Ti and T2 
denote the connected components of the forest. 

In order to define the kinematical invariants S(t'i.T2)j we introduce an arbitrary 
orientation on G. We firstly say that each external momentum pj is incoming at 
the vertex at the corresponding leg. We furthermore label each oriented edge by 
a momentum- vector g^. If the edge Ci is oriented from vertex Vj to Vk then qi is 
said to be incoming at and —qi is incoming at Vj. Momentum-conservation on 
G is reflected in our labels by the condition that the sum of all external momenta 
Pj is zero, and at each vertex, the sum of all incoming momenta is zero. By these 
conditions, except for L momenta, each of the qi can be expressed as a linear 
combination of external momenta. The kinematical invariants are defined as 



where the sign of qj is fixed by the condition that we sum over the momenta 
incoming at the component T2. Note that by momentum conservation, the S(Ti,T2) 
are functions of the external momenta. 

As an alternative to the above construction by spanning trees and forests, there 
are several ways to obtain both Symanzik polynomials from determinants of certain 
matrices [161 1131 1521, 121j . To demonstrate such a derivation, let us label each edge 
Ci by an auxiliary variable yi. Each vertex Vi is labelled by 

Zj if a leg with incoming momentum pj is attached, 
if no leg is attached. 
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Figure 1. The non-planar double-box 



For a Feynman graph with vertices Wi, Vm we consider an m x m matrix M 
whose entries are: 

_ ( Ui + J2yk for i — j, Ck attached to Vi at exactly one end, 
~ \ ~Y.yk for i ^ j, Cfc connecting Vi and Vj. 

We compute the determinant 

V(yi, .., y„, zi, Zr) = det(M) 

and consider the function 

n 

W(q;i, .., a„, zi, Zr) = V {a'^^ , .., zi, Zr) tti 

1=1 

which is a polynomial in the a— and z— variables. Note that M depends on a 
chosen ordering on the vertices but W does not. 

Let us assume that at least two legs are attached to the graph, i.e. r > 2. We 
expand W as W = W^^^ + W^^^ + ... + W'^'') where W'^'') is homogeneous of degree 
k in the z— variables. We can directly read off the first Symanzik polynomial from 
the first term in this expansion, as it satisfies 

r 

W^'^\ai, .., a„, zi, Zr) = Ug(,cxi, .., a„)^Zj. 

The massless second Symanzik polynomial 7^), g is directly obtained from VV'-^-' . By 
construction, VV(^) is homogeneous of degree 2 in the z— variables. We replace each 
product ZiZj in W*-^^ by the scalar-product of the corresponding external momentum 
vectors pi ■ pj. By momentum-conservation, J2i=iPi — Oj express each of the 
scalar products by the functions S(^Ti,T2)- result we obtain Jx),G { 161 . 

As an example let us compute the two Symanzik polynomials of the massless 
non-planar double-box, shown in figure [T] with auxiliary y— and z— variables. For 
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this graph and a chosen ordering on the vertices we have 



M = 



where Mn 
M55 = yi + 



V 

yi+y2+zi, 



Mn 

-y2 
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-y2 
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M22 


-2/3 
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M33 


~y7 





-2/5 





-yr 


M44 


-ye 











-ye 




-y4 





-y5 





-yi 


Mee 



I 



yi' 



M22 = y2 - 

= 2/4 + 2/5^ 



ye + yr + Zi, 



■y3 + z2, A/33 = ys+yb+yr, M44 

yQ, Mee = J/4 + 2/5 + ^3- Proceeding in the described way we compute 

(Zi +Z2 + Z3 + Z4)-lW(l) 

(ai + 0^2 + a3)(a4 + + + aj) + {a^ + a5){ae + ar), 

— Pi 0^2 ("1 (0^4 + as + Q!6 + a?) + (^4ae) 
-P2"2(a3(a4 + as + "6 + a?) + 050:7) 
— P3(a4a5(ai + a2 + a3 + ae + a?) + a3a4a7 + aiasae) 
— P4(a6a7(ai + a2 + a3 + a4 + a^) + a\a4a'j + a3a5a6) 
-P2)^(aia3(a4 



as + ag + ar) + aiasar + a3Q;4a6) 



-{Vi + Vzf oi2a^af, - (p2 +P3)^a2a4a7. 

It is often sufhcient to consider the Feynman integral after setting some of its 
legs on-shell, which means that the corresponding external momenta are fixed by 
setting their square to a squared particle mass. In our example we may assume 
massless particles and set = for all i — 1, 4. The corresponding Feynman 
integral was evaluated in dimensional regularization by classical polylogarithms in 
reference |59j . We will return to Symanzik polynomials of graphs with four on-shell 
legs in section m 

Let us now turn to iterated integrals. Let k be the field of either the real or 
the complex numbers and M a smooth manifold over fc. We consider a piecewise 
smooth path on M, given by a map 7 : [0, 1] — )• M, and some smooth differential 
1-forms wi, ijJn on M. The iterated integral of these 1-forms along the path 7 is 
defined by 

/ a;„...a;i = / fn{tn)dtn....fi(ti)dti, 
Jo<ti<...<t„<i 

where fi{t)dt — ^*{ijJi) is the pull-back of to [0, 1]. With the term iterated 
integral we will more generally refer to fc-linear combinations of such integrals. 

We will consider classes of iterated integrals which define the same function 
for any two homotopic paths. Such integrals are called homotopy invariant. They 
are well-defined functions of variables given by the end-point of 7. In such iterated 
integrals the differential forms and the order in which we integrate over them have 
to satisfy a property known as the integrability condition. The condition is best 
formulated on tensor products of 1-forms over some field K <Z k, which we denote 
by [a;i|...|a;m] ■ Let D denote a X-linear map from tensor products of smooth 1-forms 
on M to tensor products of all forms on A/, given by 



Z) ([a;i|...|a;„i]) = ^ [a;i|...|a;i_i|da;j|wi+i|...|w„] -I- [a;i|...|wi_i|wi A w^+il 
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A if-linear combination of tensor products ^ — X)/" o "^ii- '! ["^ii l---!"^'!]! ^i\:---.H G 
if, is called an integrable word if it satisfies the equation 

= 0. 

Let be a finite set of smooth 1-forms and let Bm{^) denote the vector space 
of integrable words of length m with 1-forms in fl. Now we return from words to 
integrals by considering the integration map on integrable words: 



[—Oii,...,iL l—Oii,...,ii ^ 

A fundamental theorem of Chen |27j states that this map is an isomorphism 
from B„i{fl) to the set of homotopy invariant iterated integrals in 1-forms in il of 
length less or equal to m, if satisfies further conditions which we do not specify 
here. 

In the following we fix X = Q and discuss two sets of 1-forms for which the 
theorem applies. For a coordinate ti on an open subset of C we firstly consider the 
set of closed 1-forms 

Ui ' ti-1' hh-r nr=i*»-i J ■ 

As a trivial consequence of dti A dti = 0, any tensor product of 1-forms in 
fl^^P is an integrable word. By applying the integration map eq. 12.31 to these 
words, we obtain the class of hyperlogarithms |46j . In particle physics it is very 
common to use sub-classes of hyperlogarithms. As an example, we may consider 
D^^^ and fix the constant t2 = —1- To physicists, the iterated integrals obtained 
from this restriction are well known as harmonic polylogarithms [54] and suffice for 
the evaluation of many Fcynman integrals. 

We want to focus on a class of functions of severable variables, obtained from 
another set of 1-forms, where now all the ti, t„ are considered to be coordinates 
in an open subset of C" : 

= 1^, where l<a<b<n 

I tl tn [[a<i<b^i ~ 1 

For this set the integrability condition is not trivial and there are words for 
which it is not satisfied. The homotopy invariant iterated integrals which we obtain 
via the integration map from the the integrable words in ^^^^ form the vector space 
B (fin) of multiple polylogarithms in n variables. We use the notation Bm i^n) for 
the vector space of such functions obtained from integrable words of length < m. 
There is an explicit map |14j to construct all integrable words in V,^^^, closely 
related to the 'symbol' in [34l IM [32] , 

The functions in were extensively studied in reference [19 . We just 

want to recall a few statements which are relevant for the following considera- 
tions. Firstly, the multiple polylogarithms of Goncharov |35) . frequently used in 
the physics literature, are contained in this class. As we want to use the elements 
of B (fin) in an iterative integration procedure, it is important for us to know their 
primitives and limits. It is proven in |19j that B (iln) is closed under taking prim- 
itives. Furthermore if we take the limits of elements of B (r2„) at i„ equal to and 
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1, we obtain Z-linear combinations of elements in B (Qn-i) , where Z denotes the 
Q-vector space of multiple zeta values. 

Now let us consider definite integrals of the form 



where / is a polynomial and l3{{gi}) £ is a multiple polylogarithm whose 

arguments are some irreducible polynomials gi. Let us call / and the gt the critical 
polynomials of the integrand. If / and the gi are linear in <„ we can evaluate the 
above integral and from the mentioned properties it is clear that the result will 
be a Z- linear combination of elements in Bm (f^n-i) . If the result can be again 
expressed by functions of the form of the above generic integrand and the critical 
polynomials are linear in then we can continue and integrate over this variable 
from to 1, and so on. 

Such an iterative procedure can be used to compute Feynman integrals. For re- 
cent examples in the physics literature, partly relying on different parametrizations, 
we refer to |30L I31L [H 1261 [5] . Aiming at such a computation one has to express 
the Feynman integral by a finite parametric integral such that the integrand can be 
written in the above form, where after each integration step, the critical polynomials 
are linear in at least one of the remaining parameters. The method was introduced 
systematically in |20j and demonstrated for certain Feynman parametric integrals 
of the type of eq. 12. 2| coming from primitive logarithmically divergent vacuum 
Feynman graphs. However, the approach is not restricted to such graphs. Refer- 
ence |23j presents a method to express Feynman integrals with UV sub-divergences 
by finite parametric integrals to which the approach may apply. The treatment of 
graphs with infrared divergences is not excluded in principle, but we are missing 
a canonical method to express IR-divergent integrals by finite ones. In principle, 
the method of sector decomposition 1 1 1 j allows us to write down the coefficients 
of a Laurent expansion for a dimensionally regularized, infrared divergent integral 
in terms of finite integrals over Feynman parameters, however, the polynomials in 
these integrals usually become very complicated. In view of the above approach 
one would ideally wish for a method, where the critical polynomials in the finite 
integrals could be obtained from the Symanzik polynomials in a rather simple way. 

For the following discussion let us assume, that in some way we have already 
been able to express a given Feynman integral by finite integrals of the type / and 
that the critical polynomials are the Symanzik polynomials of the graph. We focus 
on the criterion, that after each integration over a Feynman parameter, the new 
critical polynomials have to be linear in a next Feynman parameter. The reduction 
algorithm to be reviewed in section [3] allows us to study this criterion as it computes 
for each integration step a set in which the critical polynomials are contained. 

As a further motivation of the following discussion, let us have a glance at two 
well-known Feynman graphs in view of the mentioned criterion. For the massless 
two-loop graph of figure [2] (a) it was proven by use of the Mellin-Barnes approach 
and expansions by nested sums that each coefficient of the e-expansion is a com- 
bination of multiple zeta values |10j . Reference |20| confirmed this statement for 
this two-loop graph and several higher-loop graphs by relating them to integrals of 
the type of eq. l2.2l whose integrands satisfy the criterion. 

The case of the equal- mass two- loop sunrise graph, shown in figure [2] (b) , is 
very different. The desired coefficients in the e-expansion can be derived from the 
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m 



Figure 2. (a) Massless two- loop graph, (b) Equal- mass sunrise graph 

D = 2-dimcnsional version of the Feynman integral, 
-^^sunrisc = / daida2da3S 1 - 

with the second Symanzik polynomial 

J-G — —p^aia2a3 + inn?{aia2 + 020-3 + aia3)(ai -^02 + 03), 

playing the role of the critical polynomial. As J-q is not linear in any of the 
Feynman parameters, integral /sunrise fails the criterion. It is in perfect match 
with this simple observation, that the known result of the sunrise integral involves 
elliptic integrals f45l . In general, it is possible that the polynomials of an integrand 
fail the criterion and still the integral can be expressed by multiple polylogarithms. 
However, the criterion may provide a useful first classification and give a hint where 
to look for integrals, which exceed the class of multiple polylogarithms. We also 
want to mention reference [25], where different criteria are used to search for such 
integrals. 





3. Linear Reducibility 

Let us briefly review the polynomial reduction algorithm of [20j . Let S — 
{111 /at} be a set of polynomials in the variables ai, a„ with rational coef- 
ficients. 

(1) If there is an index 1 < ri < n such that all polynomials in S are linear in 
we can write 

fi — giOri + hi for all 1 < i < A^, 
where gi = and hi = fi\a^^=f). We define 

(3-1) ^(Vi) = {[9i)l<r<N ' {hi)l<i<N ' ih^9j - 5«/li)i<i<j<Ar} 

and furthermore we define ^(rj) to be the set of irreducible polynomials in S'^Vi)- 
S(^ri} we neglect all constants and monomials. 

(2) If there is a 1 < r2 < n such that all polynomials in 5'(ri) are linear in 
ar2, we repeat the above step, now with 5'(ri) and in the roles of S and 0;^, 
and obtain a new set of polynomials which we call S[ri]{r2)- Then, assuming that 
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starting from S the above steps can be done first for ctra and then for , we 
compute S[r2](ri) and take the intersection of both sets: 

S[ri,r2] — ^[ri](r2) ^ '5'[r2](ri) • 

Whenever we speak of intersections here and in the following, we mean the common 
zero loci, such that if a polynomial appears in two sets with a different constant 
prefactor, it nevertheless belongs to the intersection. Then we choose a next variable 
in which all polynomials of S[ri,r2] are linear and continue in the same way. At 
each iteration we apply step (1) and take the intersection 

If a set S^ri, ...,fi, ...,rk] contains a polynomial which is non-linear in a^, the set 
Sin, ...,ri, ...,rk]iri) is Undefined and omitted in the intersection. If this happens for 
all 1 < i < fc the set S[r-^,r2, ....r-k] is undefined and the algorithm stops. Unless this 
situation occurs, we obtain for each sequence of variables (an, CKr^), fc < n, a 

sequence of sets 5(^1), S[rj^,r2]i 'S'[ri, ra, r-fc] • 

Definition 3.1. We say that S is Fubini reducible (or linearly reducible) if 
there is an ordering of all n variables (a^, Q?r„) such that every polynomial in 
S[ri,r2, ...,rk] linear in Oij.^^-^ for all 1 < fc < n. 

The linear reducibility of the set of critical polynomials of an integrand as 
in eq. 12.41 is a criterion for the integral to be computable by the above approach. 
The criterion is sufficient but not neccessary. The sets S^^.-^^ .... r-fc] contain the critical 
polynomials of the integrand after the first k integrations, but might as well contain 
spurious polynomials which drop out in the integration procedure. A more refined 
reduction algorithm presented in reference |21j omits such cases, but not necessarily 
all of them. Furthermore the occurrence of a quadratic polynomial does not always 
forbid us to continue with the computation. 

By applying the above algorithm to first Symanzik polynomials, it was shown 
in '20] that several vacuum Feynman integrals can be computed and evaluate to 
combinations of multiple zeta values. Moreover the same is true for coefficients of 
a dimensional series expansion of Feynman integrals with two legs, obtained from 
cutting one edge in one of these vacuum graphs, as in the case of figure [2] (a) . The 
linear reducibility of first Symanzik polynomials is extensively studied in terms 
of Dodgson polynomials in [21 . It is shown that the first five iterations of the 
reduction succeed for any first Symanzik polynomial. Moreover, a first Symanzik 
polynomial is reducible, if its graph has vertex width less or equal three. (This is 
the class of graphs which decompose into two connected components after removing 
three vertices or less.) These results explain, why one has to go up to complicated 
graphs at high loop orders to find first examples for vacuum-type Feynman integrals 
which exceed the set of multiple zeta values |24l I22j . 

In the following, we want to consider the above algorithm applied to both 
Symanzik polynomials. In order to include which may depend on particle masses 
and external momenta, we slightly extend the above formulation of step (1), allow- 
ing fi to be polynomials whose coefficients are rational numbers or algebraic func- 
tions of additional parameters si, Sm- The rest of the algorithm is not affected 
by this change. 

It will be useful to consider polynomial reduction in coordinates for products 
of P^. Let P{ai, an) be a polynomial in n Feynman parameters and consider 
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an associated function of new variables xi, a;„, yi, ?/„ defined by 

where di{P) denotes the maximal degree of P in for all 1 < i < n. We obtain 
our original polynomial back by setting 

(3.2) -P(ai, a„, 1, 1) = P{ai, a„). 

The function P is by definition a polynomial in xi, x„, yi, y„ and it is linear 
in Xk if and only if P is linear in a^. 

Let us consider a set of polynomials S — {Pi, Pn} which are all linear in 
afc. We can write 

d 

Pi = Oik^ Pi + Pi\ak=0- 

oak 

Consider the corresponding set of polynomials S — {Pi, Pn}, obtained from S 
by changing to the projective coordinates. Each of these polynomials is linear in 
Xk and we can write 

Pi = XkPi\xk = l,yk=0 + ykPi\xk=0,yk = l- 

Let us use the convention that in a Fubini reduction with respect to the x-variables, 
we neglect the above prefactor y^ in the sense that in eq. 13.11 the terms are given 
by 5i = Pi\xk=i.yk=o and hi = Pi\xk=o.yk=i: which does not affect the linear re- 
ducibility. Step (1) of the reduction algorithm applied to S with respect to Xk then 
gives the set S(^k) consisting of the irreducible factors of 

^(k) — {(-PiUfc = l,yfc=o)i<j<7v ' (-^*Ufc=0,iyfc = l)i<i<Ar J 

(3.3) [Pi\xk=a,yk = l ' Pj\xk=l,yk=0 ~ Pi\x^ = \,yk='A ' Ulc=0, i/fc=l) | ■ 

Considering a Fubini reduction in x-variables instead of a-variables, we con- 
vince ourselves that the factorizations into irreducible polynomials are not affected 
by the change of variables. Indeed if P factorizes as 

P = fl-f2 

into two polynomials /i , /2 then P factorizes as 



/l-/2 





Vn 



because of di{fi) + di{f2) = di{P) for all 1 < i < n. As a consequence we obtain: 

Lemma 3.2. P is linearly reducible with respect to the variables ai, a„ if 
and only if P is linearly reducible with respect to the variables Xi, a;„. 

Now let us choose some coordinate xi, 1 < I < n, and define two new sets of 
polynomials by restrictions x; = 0, y; = 1 and xi = 1, yi = respectively: 
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To these sets we apply one reduction step with respect to Xk- At first let us 
assume I ^ k. Step (1) of the algorithm gives 

'^(i)"^'' ~ irreducible factors of la;, =o,3/i=i) 
^(k) '^^ = irreducible factors of la;, =i_j^,=o- 
If a polynomial P factorizcs as P ~ fi ■ f2 then furthermore 

P\xi=0,vi = l = /l|a;,=0,i/, = l • /2U,=0,i/, = l, 
P\xi = l,yi=0 = fl\xi=l,yi=0 ■ f2\xi=l,yi=0, 

and the maximal degrees with respect to another variable Xi satisfy di{fj) > 
di {fj\xi=o,y,=i) and di{fj) > di {fj\xi=i.y,=o) for j = 1, 2 and for all i = 1, n. 
This means that for each polynomial / G S'(fe) there is at most one polynomial 

/' G S'^'i^;-I^''^\ and di{f) > di{f') for all i = \, n. The same is true for S^l^^^^'^\ 

Lastly we observe that in the case oi I — k, the irreducible factors of 5*^'' ^' 
and S^''^'*^' are already contained in S(k). This proves the following lemma: 

Lemma 3.3. If S is linearly reducible with the ordering (a;^, x^) then for 
any 1 < I < n the sets 5|a;,=o, 3/1=1 and 5|a;,=i,3/,=o are linearly reducible with 

(^ri 5 "^li "^rn) ■ 

In combination with the previous lemma we obtain: 

Lemma 3.4. Let S = {Pi, Pn} be a set of polynomials which is linearly 
reducible with (a^, a^^) and whose members are linear in ai. Then the sets 

~ {§^' '^'^^ ~ {Pi\ai=o, PN\ai=o} drc linearly reducible with 

(*^ri 5 •••■) *^rn)* 

4. Towards a Classification by Critical Minors 

Let G be a graph and Eg its set of edges. For e € G we denote by G\e the 
graph obtained from G by deletion of e. Furthermore we write G//e for the graph 
obtained from G by contraction of e. This is the graph where the end-points of e are 
identified and then e is removed. For any distinct edges ei , 62 € Eg the operations 
of deleting (or contracting) ei and deleting (or contracting) 62 commute. Therefore 
we can more generally write 7 = G\D/ /C with distinct D, C C Eg, for the unique 
graph obtained from G by deleting all edges in D and contracting all edges in C. 
Any such 7 is called a minor of G. 

If G is connected and there is an edge e such that G\e is disconnected then e 
is called a bridge. When speaking of Feynman graphs, we may ignore disconnected 
graphs, so we introduce the convention that the corresponding Symanzik polynomial 
are zero: 

i^G\e ~ 0, Ta.G\e = if 6 is a bridge. 

Let us furthermore call G Eg a tadpole if it is attached to the same vertex at both 
ends. In this case we have a factorization in the corresponding Feynman parameter: 
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For any e which is not a tadpole we have the well-known deletion/contraction 
identities: 

Ug = UG\e'^e +Uc/ /e^ 
J^O,G = J^Q,G\eae + J^O,G//e- 

Now we consider a set of arbitrary particle masses mi, m„ distributed over the 
edges of G, with the restriction that at least one edge is massless, mk = 0. Then 
the second Symanzik polynomial 



j^G — ^o.g+ug ■ 



is linear in at and the above relations are extended by 

Fg = J^GVefcttfe if efc is a tadpole, 
Tg = J^G\ek'^k + J^G//ek if is not a tadpole. 
Equivalently, if Cfe is a massless edge of G and not a tadpole then 
d 

(4-1) -^^Ug ^UG\eu,UG\ak=0=l^G//eu^ 

d 

(4-2) 'da'k'^'^ ^ -^G\efc, -^cUfc^O = ^G//ek- 

In the following let us call a Feynman graph G linearly reducible, if {Ug, J'g} is 
linearly reducible. We arrive at the main statement of these notes: 

Theorem 4.1. If G is a linearly reducible Feynman graph then any minor of 
G is linearly reducible as well. 

Proof. Let G be an arbitrary Feynman graph, ej any of its edges and G = 
G\m j=o the graph obtained from G by setting the mass mj associated to Cj equal 
to zero. Assume that G is linearly reducible. Then G is linearly reducible for any 
value of mj and therefore G is linearly reducible as well. As Cj is massless in G we 
can apply equations 14.11 and 14.21 and obtain, if Cj is not a tadpole, 

^Gf/e, = {^G//e,, ^Gljej] = {U^\a,=n, J^g.U,=o} • 

If Cj is a tadpole we have 

By lemma [33] it follows from the linear reducibility of G that G\ej and G//ej 
are linearly reducible. G and G have the same minors with respect to deleting or 
contracting Cj : 

G\e, = G\e„ G//e, = G//e, 

and therefore SG\ej = Sq-^^. and SG//ej ~ ^G//e - induction, this proves the 
theorem. □ 
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Figure 3. K4 with four legs 



A set of graphs Q is called minor closed if for all G S 5 every minor of G 
belongs to Q as well. Let V. be any set of graphs and let Q-^ be the set of all graphs 
which do not have a minor in the set T-L. Then the set Q■^ is minor closed, and 
the graphs in H are called forbidden minors of Q-u. A theorem of Robertson and 
Seymour [55j states that any minor closed set of graphs can be defined in such a 
way by a finite set of forbidden minors. 

A graph H is called a critical minor of a set of graphs Q if the minors of H 
belong to Q but H does not. By removing a graph from H which is a minor of 
another graph in T-L, we do not change the set G-u. Therefore, to define a minor 
closed set Q-h, it is sufficient to let T-L consist only of critical minors. A well-known 
example for a characterization by critical minors is given by the set of planar graphs. 
Due to a theorem of Wagner 61 , the set of planar graphs is the set with forbidden 
minors {^^3^3, K5}. 

In general we have to be careful when adapting notions from pure graph the- 
ory to the study of Feynman graphs, which are equipped with labels and special 
properties. However, theorem l4.ll suggests to attempt an analogous characteriza- 
tion of the set of linearly reducible Feynman graphs by critical minors, i. e. by a 
set of Feynman graphs which are not linearly reducible but have linearly reducible 
minors. We conclude these notes with a first step towards such a characterization 
for Feynman graphs with a dependence on external momenta. Let A be the set of 
massless Feynman graphs with four on-shell legs, attached to four distinct vertices. 
The on-shell condition restricts the four external momenta to satisfy pf = such 
that the dependence on the external momenta can be expressed by two Mandelstam 
variables. By use of a Maple-implementation of the Fubini algorithm, we find lin- 
ear reductions for all two-loop graphs of this class. At three loops there are several 
graphs in this class for which our program fails to find a linear reduction. Most of 
these graphs have the graph of figure[3]as a minor. This is the complete four-vertex 
graph where we attached an on-shell leg at each vertex. From direct observa- 
tion of the two Symanzik polynomials one can see that the graph is not linearly 
reducible. As on the other hand all its minors are linearly reducible, the graph 
plays the role of a critical minor for the set of all linearly reducible graphs. Our 
case study suggests that only a few further critical minors are needed to distinguish 
all linearly reducible graphs of A at three loops. In this way, a large class of graphs 
can be separated into linearly reducible and irreducible members by knowing only 
a small number of critical minors. 
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5. Conclusions 

In this talk we reviewed the criterion of linear reducibility of Synianzik polyno- 
mials, which can be used to decide whether a corresponding Feynman integral can 
be computed by iteratively introducing hyperlogarithms or multiple poly logarithms. 
Brown showed in [21j that with respect to the first Symanzik polynomial, the set 
of linearly reducible graphs is closed under taking minors and that therefore a char- 
acterization of this set by forbidden, critical minors is possible. We extended this 
line of argument to the case of graphs with masses and kinematical invariants, by 
showing, that minor closedness is true for both Symanzik polynomials. We exhibit 
a first critical minor with a non-trivial dependence on kinematical invariants. 

We expect a classification with respect to the criterion of linear reducibility 
to be useful for the more difficult question of which Feynman integrals evaluate to 
multiple poly logarithms. Due to the simplicity of its derivation and the property 
of minor closedness, such a systematic classification of a large class of Feynman 
graphs appears feasible. 
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